The challenge of unravelling magnetic properties in LaFeAsO 
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First principles calculations of magnetic and, to a lesser extent, electronic properties of the novel 
LaFeAsO-based superconductors show substantial apparent controversy, as opposed to most weakly 
or strongly correlated materials. Not only do different reports disagree about quantitative values, 
there is also a schism in terms of interpreting the basic physics of the magnetic interactions in this 
system. In this paper, we present a systematic analysis using four different first principles methods 
and show that while there is an unusual sensitivity to computational details, well-converged full- 
potential all-electron results are fully consistent among themselves. What makes results so sensitive 
and the system so different from simple local magnetic moments interacting via basic superexchange 
mechanisms is the itinerant character of the calculated magnetic ground state, where very soft 
magnetic moments and long-range interactions are characterized by a particular structure in the 
reciprocal (as opposed to real) space. Therefore, unravelling the magnetic interactions in their full 
richness remains a challenging, but utterly important task. 

PACS numbers: 



X 



The disovery of high temperature superconductiv- 
ity in LaFeAs0i_2,F 2 . by Kamihara and co-workersi 
rising to critical temperatures T c over 50K with 
rare earth substitution^ has resulted in a great deal 
of exp er iment a& 3 -^£& 7 -&£ i 1 ' 1 1 ' 1 2 and theoretical 

ac ^ v ^„ 13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29 These 

studies have shown that this set of Fe-based super- 
conductors displays a remarkably rich set of physical 
properties. Besides the high critical temperatures, which 
include the highest known values of T c outside the 
cuprates, there is evidence for several types of Fermi sur- 
face nestin g 15 i 19 which raises the possibility of typically 
metallic collective excitations, such as itinerant magneti- 
zation waves. At least three different competing types of 
magnetic fluctuations have been predicted 9 ^ 13 ' 14 ! 15 ! 16 ? 17 
and an ordered spin density wave was first predicted 9 -^ 
and subsequently observed experimentally in the 
non-superconducting undoped parent, LaFeAsQ 10 ^ 1 ' 12 . 
Currently, most researchers agree that superconductivity 
in this compound is unconventional and likely related to 
magnctic m 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24,30 - 31 



ferent DFT calculations disagree among themselves and 
that, beyond a general initial consensus on the common- 
alities of the different materials, electronic structure cal- 
culations differ in details of the ground state and in the 
band structure near the Fermi level. The reason for such 
inconsistencies is not just computational inaccuracy, as 
it may seem, but rather has a physical origin. Namely, 
as pointed out by a number of authors and discussed 
in detail below, magnetism in this compound is very 
itinerant, making the calculated magnetic energies and 
moments extremely sensitive to the approximation used 
and to tiny details of the crystal structure. This situa- 
tion is relatively unusual, particularly in comparison to 
compounds with localized magnetic moments, which are 
normally rather robust. Thus, there may be a tendency 
to ascribe reported differences to an inherent inaccuracy 
of electronic structure calculations, rather than to the 
peculiar dependencies of the compound itself. Here we 
outline those dependencies and provide a reference set of 
calculations in an attempt to establish a clear theoretical 
picture within DFT. 



The uncommon richness of the electronic struc- 
ture of this compound has already led to a situa- 
tion in which different theoretical groups report den- 
sity functional theory (DFT) band structure calcu- 
lations that emphasize proximity to different mag- 
netic states: weak ferromagnetis m 13 i 14 , checkerboard 
antifcrromagnctis m 13 ! 15 ! 17 , an antiferromagnetic stripe 
phase 9 -^. For the general reader, it may look as if dif- 



Therc is a kernel of truth underlying the perception 
that DFT results are disparate: for soft itinerant mag- 
nets one has to exercise extra caution in the calculational 
parameters in order to get consistent and reliable re- 
sults. In particular, methods that imply restriction on 
the shape of the charge density (such as Atomic Sphere 
Aproximation) and methods that involve pseudization 
of the crystal potential can be trusted only after they 
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have been tested against all-electron, full-potential cal- 
culations. Furthermore, in this particular compound the 
magnetic properties show an unusual sensitivity to the 
internal As height and are also sensitive, unusual for lo- 
calized moments but typical for itinerant magnets, to the 
exchange correlation potential. 

Before investigating the magnetic properties and their 
dependencies in detail, it is important to note that this 
system is in proximity to a quantum critical point that 
results in strong spin fluctuations and to establish how 
such a system is best dealt with theoretically. It is well 
known (Refs. 1 3 2j j 3 3l and references therein) that DFT 
calculations, which are mean-field by nature, underes- 
timate the effect of spin fluctuations that generally de- 
press or suppress long-range magnetic order. It is also 
understood that gradient corrections (GGA) to the local 
density approximation generally increase the tendency 
to magnetism. As a rule of thumb, in systems with 
strong local Coulomb correlations, such as cuprates or 
3d metal oxides, GGA better describes magnetic prop- 
erties (an even further improved description is provided 
by the LDA+U method), while in itinerant magnets or 
near-magnets, LDA is closer to experiment in terms of 
magnetism. The classical examples are Pd metal and 
Si^RuO^ both are ferromagnetic in GGA, but not (in 
agreement with the experiment) in LDA. On the other 
hand, GGA is usually better at predicting the crystal 
structure of transition metal compounds. Finally, given 
the same crystal structure, the bands themselves arc 
practically indistinguishable whether calculated in LDA 
or GGA. With these facts in mind, previous experience 
in first principles DFT calculations suggests that the best 
approach to the LaFeAsO family may be to optimize 
the structure (if necessary) using GGA and to calculate 
magnetic properties using LDA (while LDA consistently 
overestimates the tendency to magnetism compared to 
the experiment, it always fares better than GGA). Ei- 
ther functional can be used to analyze bands and Fermi 
surfaces. If a quantitative value of the magnetic stabiliza- 
tion energy is less important than an accurate analysis 
of a trend (say, with doping), GGA may be chosen for 
the reason that it amplifies the tendency to magnetism 
and results in larger numbers that are easier to com- 
pare with each other. What is important though is that 
as long one uses full-potential, all-electron methods (or 
well tested pseudopotentials) the results agree satisfacto- 
rily with each other, within the same density functional 
(LDA or GGA). It is worth noting that one can envision 
a situation where proper account for the energy associ- 
ated with magnetic fluctuations may be more important 
than one-electron excitation spectrum (even when ad- 
dressing the mechanism of superconductivity). In such 
a case, spin-polarized GGA may give more insight than 
LDA (despite yielding an obviously exaggerated magnetic 
ground state) ; there is little previous experience to guide 
band theory procedure in this case. 

Let us now discuss the reported results of magnetic 
calculations in more detail. Three different long range 



TABLE I: Magnetic stabilization energies (vs. non-magnetic) 
and magnetic moments for zero doping and hole (x = 0.1) 
doping in the experimental and optimized structures. Unfilled 
boxes correspond to unstable (vs. metastable) configurations. 
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magnetic orders in the Fe plane have been so far consid- 
ered: a ferromagnetic ordering that retains the crystal 
symmetry P4/nmm (#129) with two formula units per 
cell, a checkerboard antiferromagnetism (P4m2, #115, 
two formulas), and an antiferromagnetism with alternat- 
ing stripes along 100 (Pccm, #49, four formula units). 
Singh and Dr>H, following the prescription above, pre- 
sented LDA calculations for the undoped and 10% doped 
materials. They found, using an in-house full-potential 
LAPW code, that the undoped material is barely unsta- 
ble against weak ferromagnetism (M < 0.1/is/Fe), but 
stable against checkerboard antiferromagnetism. Simi- 
larly, it was found that in LDA the ferromagnetic insta- 
bility rapidly disappears with doping and at x =0.1 no 
longer exists. Mazin et al^ later discovered that within 
LDA the doped material is unstable against the stripe 
AFM ordering. It was shown subsequently shown£ that 
this is the universal ground state within DFT at all dop- 
ings including zero. 

Soon after, Cao et al, using pseudopotential meth- 
ods, reported GGA results for the checkerboard AFM 
ordering in the absence of doping. They used two 
different techiqucs, VASP+PAW^ and PW-SCF ultra- 
soft pseudopotential^. Neither of these methods is all- 
electron and therefore depends on the selected pseudopo- 
tential (especially for Fe). Not surprisingly, while qual- 
itative conclusions agree with all-electron calculations, 
quantitatively they differ both between themselves and 
with the all-electron "exact" result. Their VASP calcu- 
lations gave a very large magnetic stabilization energy of 
84 mcV/Fe, while PW SCF produced 14 meV/Fe. 

Similarly, Dong et a& reported pseudopotential calcu- 
lations for the stripe ordering proposed in Ref.— , also in 
GGA, and found a stabilization energy of 40 mcV/Fe and 
magnetic moment of 1.5 /is (Dong et al do not indicate 
what level of doping they were using, but their plots are 
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FIG. 1: (color online) Fixed spin moment calculations for 
LaFeAsOi-^Fa. 



consistent with x=0.1). 

An extensive study of both types of AFM ordering was 
reported by Yildirim^ who performed both pseudopo- 
tential (which he discarded as unreliable) and all-electron 
calculations using LDA. Using fixed spin moment calcu- 
lations, he found that only the stripe order is stable with 
respect to a non-magnetic state. Ma et aU&, using GGA 
and a pseudopotential approach, also found the stripe or- 
der to be the ground state, but also found a metastable 
checkerboard AFM state. 

Finally, Yin et a&- performed GGA calculations and 
observed that (a) the magnetic moment depends on the 
As position anomalously strongly and (b) the As posi- 
tions, optimized in GGA (as it is generally considered to 
be the most accurate for the structural properties), differ 
from the experimentally reported by as much as 0.1 A, 
an exceptionally large discrepancy. 

With this in mind, we have performed calculations us- 
ing a variety of codes: two entirely different all electron 
full-potential codes (WIEN2k^ and FPLQ2I), as well as 
two different pseudopotential codes (VASP— and PW 
SCF—). The two all-electron codes gave nearly iden- 
tical results, while the pseudopotential results differed 
somewhat, though not drastically, both from all-electron 
and from each other to some degree. We pursued the 
"standard" route discussed in the earlier section (using 
GGA for optimization and then recalculating the mag- 
netic properties using LDA). Optimization is performed 
in the paramagnetic phase, since experimentally these 
materials are either paramagnetic or have a drastically 
suppressed magnetic moment of 0.15-0.35 [Lb- The ex- 
perimental lattice parameters, a =4.035 and c= 8.741 
A were used. We also compare our results to those in 
LDA (in the same stucture) and to both LDA and GGA 
calculations in the experimental structure. 

Optimization leads to As and La positions that vary 
little with doping: from z^ a = 0.1450, za s = 0.6380 at 



x = to z La = 0.1490, z A s = 0.6340 at x = 0.2. Reported 
experimental data are, respectively, 0.1415 and 0.6513 at 
zero doping, and practically the same at x = 0.1 (ref.). 
As pointed out by Yin22, the discrepancy in the As posi- 
tion amounts to 0.1 A, a nearly unheard-of an error for 
GGA calculations. In Table [J and Fig. [T] we show how 
the magnetic stabilization energy varies as a function of 
doping for the optimized and experimental structures. 
One sees that all GGA results are substantially more 
magnetic (expectedly) as are the calculations in the ex- 
perimental structure (unexpectedly). A consistent set of 
numbers for calculations using all-electron WIEN2k are 
given, with comparison numbers from other codes also 
included. The two all-electron methods differ by about 
10% in terms of magnetic energy; pseudopotential calcu- 
lations are somewhat more off, especially the PW SCF. 
Importantly, all of the them predict the same energy se- 
quence (this is however not the case if a pseudopotential 
without p states in the valence is chosen for Fe2£). Fi- 
nally, we compare the results with published calculations 
and find the same trend: all electron calculations agree 
reasonably well among themselves, while pseudopoten- 
tial ones, being sensitive to the choice of pseudopotential, 
scatter more. 

Interestingly, apart from the structural parameters 
themselves LDA results obtained using the GGA- 
optimized structure seem to correspond more closely to 
what is physically observed than those obtained using 
the experimental structure. General experience indicates 
that magnetism in metals near quantum critical points 
can be entirely suppressed by fluctuations if the LDA 
magnetic energy is on the order of 10-15 meV per atom, 
but it is rather hard for spin fluctuations to entirely de- 
stroy magnetism that is stabilized by any substantially 
larger energy. The stabilization energy for the stripe anti- 
ferromagnetism in the experimental structure is 75 meV 
(x = 0.1), far too much to expect full suppression of 
magnetic order from spin fluctuations. In the theoreti- 
cal structure, on the other hand, the same energy is 9 
meV, which is in precisely the range of energies that ac- 
cumulated knowledge of itinerant magnets indicates can 
be overcome by spin fluctuations. 

The Fermi surfaces corresponding to the GGA- 
optimized structure in comparison to those of the exper- 
imental structure further support use of the optimized 
structure. An empirical observation can be made that 
the border between the SDW antiferromagnetism and 
superconductivity in LaFcAsOi-^F^ is roughly at the 
same concentration (x ~ 0.03)^ where the only substan- 
tially 3D piece of the Fermi surface disappears (x ~ 0.05). 
At the same doping level, magnetic ordering disappears, 
very much in line with standard intuition that long range 
order in a 2D system is easily destroyed by fluctuations. 
Superconductivity appears at the same time, suggesting, 
rather naturally, that at smaller x it is simply suppressed 
by the magnetic order. The nice correlation between di- 
mensionality, magnetism, and superconductivity is en- 
tirely destroyed in calculations based on the experimcn- 
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FIG. 2: (color online) A comparison of the band structure 
derived using experimental coordinates to that using opti- 
mized coordinates. The three-dimensional band that crosses 
the Fermi energy in the optimized structure is extraordinarily 
sensitive to As height, dropping below Ef when it is raised to 
the experimental position. It is replaced by a two-dimensional 
band. 



tal structure. The 3D band that crosses Ep and forms 
a small pocket around the Z point of the BZ (see Fig. 
[5]) in the optimized structure is far below Ep in the ex- 
perimental structure, replaced by a strongly 2D band. 
In other words, expanding the Fe-As distance exchanges 
the positions of a 3D and a 2D band, causing the latter 
to cross Ep. The system is then two-dimensional for all 
values of x and doping is expected to have little effect on 
the strength of the spin fluctuations. Of course, at some 
doping this 2D band will also fill (similar to the 3D band 
of the optimized structure), but no dimensional crossover 
will occur coincident with the disappearance of the Fermi 
surface. 

Although the experimental trends can seemingly be 
better explained by using the GGA-optimized coordi- 
nates followed by LDA magnetic calculations than by 
using GGA magnetic calculations or experimental coordi- 
nates, as outlined above, the exception of the coordinates 
themselves is important. Surprisingly, applying spin po- 
larization and the GGA exchange-correlation potential 
during optimization leads to rather good agreement with 
experiment for internal positions as well as lattice con- 
stants (the latter are already reasonably reproduced in 
paramagnetic GGA calculations). The Fe-As height, to 
which the FS and magnetism are unusually sensitive, is 
obtained as 2.38 A, acceptably close to the experimen- 
tally observed height of 2.41 A. Thus, there is a contra- 
diction between good agreement in structural parameters 
and wildly overestimated magnetism (or conversely, good 
magnetic energies/moments and badly underestimated 
bond-lengths) that remains one of the more confusing 
aspects of this compound. In order to fully understand 
how magnetism manifests itself in this family of materi- 



als, which is crucial as a beginning step to understand 
the relationship between magnetism and superconduc- 
tiviy, the connection between the physical and magnetic 
structures must be disentangled. This remains an open 
and intriguing problem. 

We now address the origin of the antifcrromagnctic 
interactions. As pointed out in RefJ^, all three mag- 
netic instabilities have entirely different physical origins: 
ferromagnetism is of the Stoner type, checkerboard anti- 
ferromagnetism results from the combined effort of con- 
ventional superexchange and of nesting of two electron 
pockets^, while the stripe ordering appears because of 
nesting between the hole and electron pockets, and is ad- 
ditionally supported^ by the next-n.n. superexchange. 
It has been assumed by some author o 28 i 30 ' 31 i 40 that su- 
perexchange is the leading cause of both AFM insta- 
bilities, and, by implication, that a two-shell Heisen- 
berg model provides a meaningful description of mag- 
netism. Other o 9 ' 15 i 19 emphasized the long-range charac- 
ter of magnetic interaction, coming from interband spin 
susceptibility The physical difference between the two 
descriptions is enormous: the former implies AFM in 
real space, that is, each two individual ions that are neig- 
bors or next neighbors to each other interact antiferro- 
magnetically, but there is no direct interaction at longer 
distances. The latter interaction is local in the recipro- 
cal space, that is, it corresponds to condensation of a 
spin-density wave (even though only a finite amplitude 
SDW is stable), and requires ordering on the scale of 
at least several unit cells to exhibit antifcrromagnctism. 
Direct evidence in favor of the long-range SDW scenario 
comes from the calculations of Yin et. aS^\ they have 
computed nn and nnn exchange constant in the ordered 
stripe phase with respect to the small deviations from the 
AFM ordering. They found that the exchange constants 
have opposite signs along the CDW wave vector and in 
the perpendicular direction. Given that the underlying 
structure is tetragonal and the x/y asymmetry arises en- 
tirely due to the magnetic ordering itself, this proves be- 
yond any doubt that the Heisenberg Hamiltonian is not 
applicable for this system. 

We have also undertaken a more conventional test that 
also confirms the SDW character of the stripe antifcr- 
romagnctism. First, one can try to map the magnetic 
energy of the two stable AFM structures onto the n.n. 
Heisenberg model. This is a routine procedure for lo- 
calized systems. Here, however, it has the problem that 
the moments are so soft that not all desired magnetic 
patterns can be realized, most notably the FM struc- 
ture, and different patterns converge to different mag- 
netic moment amplitudes (extreme softness of the mag- 
netic moment magnitude already implies that the appli- 
cation of the Heisenberg Hamiltonian is suspect). This 
problem was realized by Yildirim who added an exter- 
nal field (when needed, staggered) to converge all three 
patterns he considered at the same moment magnitudes. 
The hope is that the result is not too sensitive to the 
shape of the staggered field (which is true for localized 
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moments, but not necessarily for itinerant ones). Ma, 
on the other hand, allowed the moments to vary freely 
within each of the three antiferromagnetic orderings, and 
found magnetic moments ranging from 2.2 - 2.6 /is, and 
then adjusted the ferromagnetic structure to have ap- 
proximately the same moment. 

In the checkerboard structure, the Heiscnbcrg 
(or Ising, in this case) exchange energy per site, 
Y,i>j JijMiMj, is E cb = -2(J n „ - J„„„)M 2 , and in the 
stripe phase E st = —2J nnn M 2 (with respect to the non- 
magnetic state). It is also possible^ to stabilize an in- 
termediate magnetic structure, where the sites with the 
coordinates (2n, 2m) have spins up, (2n + 1, 2m + 1) 
spins down, and (2n + 1,2) and (2n, 2m + 1) are not 
magnetic. In that case the magnetic energy per site is 
half that of the stripe phase, — J nrm M 2 . Finally, the 
ferromagnetic state energy in the Heiscnbcrg model is 
+2( J nn + J nnn )M 2 . Ma et al^ give for the ferromagnetic, 
checkerboard, stripe and the intermediate structures the 
energies of +91, —217, —109 and —65 meV/Fe, respec- 
tively. It is obvious that these four numbers are nowhere 
close to be mappable onto this model, contrary to the 
claim in Ref4£. It is also worth noting that the very fact 
that the checkerboard structure is stable in GGA implies 
Jnn > Jnnm contrary to the popular belie f 30 i 31 that the 
nnn superexchangc is stronger here. 

One may notice that the energies above are inconsis- 
tent with any published all-electron calculation, so one 
may agree with Yildirim 2 ^ that in this case pseudopo- 
tential calculations cannot be trusted (although our own 
experience indicates that carefully and properly selected 
pscudopotcntials are still reasonable). To test that, we 
have performed similar calculations using the all-electron 
WIEN code, also using the optimized structure, x = 
and GGA (which makes sense in this particular case in 
order to artificially emphasize the magnetism) . We found 
the four energies in question to be +166, —18, —81 meV 
and —36 meV (and the moments ranging from 1.5 to 
1.8 /j,b)- These numbers are reliable DFT results, with 
no uncontrollable approximations such as pscudization 
of the potentials, yet they are just as incompatible with 
the two- neighbor Heiscnbcrg model as are Ma et al's re- 
sults, in agreement with the perturbative calculations of 
Yin et. al. 

We conclude this section by restating the main points 
relevant to magnetic calculations in this system. In both 
GGA and LDA the undoped system is on the verge 
(within a fraction of a meV) of an instability against 
an itinerant ferromagnetic state with a very small mag- 
netic moment. The AFM checkerboard structure, remi- 
niscent of the magnetic ordering of Cu planes in cupratcs, 
is quite stable, albeit not the ground state, in both LDA 
and GGA, if one uses experimental As positions, in ci- 
ther the doped or undoped compounds. In the optimized 
structure the checkerboard magnetism is stable only in 
GGA, and by a rather small energy. The striped struc- 
ture is the DFT ground state for all relevant dopings, 
independent of the As position; in GGA it is very star 




FIG. 3: (color online) GGA Fermi surfaces of LaFeAsOi-ajFz, 
for a) x—0 with optimized coordinates, b) x=0.1 with opti- 
mized coordinates, c) £=0 with experimental coordinates, and 
d) x—0.1 with experimental coordinates. Note that the band 
structure is three dimensional at low doping in the optimized 
structure, but becomes two-dimensional with doping, while in 
the experimental structure it is two-dimensional and insensi- 
tive to doping 

ble, while in the LDA, and if one uses optimized po- 
sitions, it rapidly becomes a borderline instability that 
can be easily destroyed by fluctuations, in agreement 
with the experiment. Magnetism cannot be described, 
not even qualitatively, using two superexchange coupling 
constants, not even if one adds a Stoncr energy reflecting 
the softness of individual spins. Only LDA calculations, 
and only in the optimized structure, provide a quantita- 
tively realistic picture of the magnetic properties of this 
system in the entire doping range. It is worth reiterat- 
ing that in borderline magnets, of which this system is 
clearly an example, not only is LDA closer to experiment 
than GGA, but in real materials fluctuations suppress 
magnetism even further, so that in comparison to GGA, 
LDA is a step in the right direction, but does not go the 
whole way. 

In the second part of this work, we concentrate on a 
detail of the Fermi surface that is crucially important for 
some of the proposed superconductivity theories, as well 
as for some of the magnetic orderings: the eccentricity 
of the overlapping ellipses around the M point of the 
BZ. The importance of the eccentricity is obvious from 
the fact that if the electron pockets were ideal circles 
(though they do not have to be, by symmetry), just as 
the hole pockets, they would nest ideally with each other, 
creating a sharp nesting peak at Q = 7r, 7t and strongly 
enhancing the tendency to checkerboard AFM ordering. 
At the particular doping where the number of holes is ex- 
actly equal to the number of electrons, the hole and the 
electron pockets would also nest ideally, making this spe- 
cific concentration very special from the electronic point 
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TABLE II: Calculated eccentricities, e(F) and e(Z) of the 
electron Fermi surfaces as a function of doping level with dif- 
ferent methods. Here e(r) and e(Z) are in the k z =0 and 
k z =0.5 planes, respectively. Note that at low doping there is 
an additional 3D hole sheet. 
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of view. In addition, such sharp structure in the mag- 
netic susceptibility is favorable for exotic order param- 
eter distributions with nodes, set not by symmetry but 
by resonance effects between the peaks in the suscepti- 
bility and the Fermi surface geometry-^. At least one 
of the proposed theorie s] 20 ' 21 relies on the condition that 
two rotated electron pockets coincide with an accuracy 
in energy not worse than the superconducting gap. 

Although the band structure and fermiology are much 
less sensitive to calculational method than the stabiliza- 
tion energies are, there is still some variation between 
methodologies and substantial differences occur between 
published reports (in the latter case, largely due to vari- 
ance among the structures used by different authors). 

We have summarized our own results, again using two 
all electron full-potential codes and two pseudopotential 
codes in Table UH where we additionally include eccen- 
tricities extracted from other reported calculations. The 
eccentricity in the T—X—M (k z = 0) plane is listed, and, 
where available, in the Z — R — A plane (k z = 0.5). The 
latter is universally larger, due to the fact that the longer 
axis of the ellipse is formed by the only band around the 
M point that is (somewhat) z-dispersive. 

Again we see excellent (nearly perfect) agreement be- 
tween our two all electron results with more deviation in 
the pseudopotential ones. In fact, all pseudotential cal- 
culations appear to underestimate the eccentricity, some 
by a large amount. Ref. [2l| has reported zero eccentric- 
ity at a doping of ~ 10% in the T plane, far less than 
not only both all-electron calculations, but also other 
pseudopotential results. Surprisingly, in the Z plane, the 
eccentricity very closely matches our pseudopotential re- 
sults, indicating a highly dispersive band absent from 
our calculations. This is in some sense unfortunate, be- 
cause a very interesting model of interband triplet s-wave 
pairing 2 ^ is only viable if the energy difference between 
the two bands at the M point is smaller near the point 
where they cross the Fermi level, than the superconduct- 



ing gap. Our all-electron results unequivocally indicate 
that this difference is at least 100 mcV in the k z = 
plane and at least 180 mcV in the k z — ir/c plane. 

Several things should be noted when comparing these 
calculations. The first is that doping is accounted for dif- 
ferently in our pseudopotential calculations than in our 
all-electron calculations. In the former, the desired num- 
ber of electrons are added and compensated for with a 
uniform positive background. In the latter, electron dop- 
ing is obtained using a fictitious O atom with increased 
core charge and electrons and hole doping is obtained us- 
ing a fictitious La atoms with decreased core charge and 
electrons. Thus, one should not expect the two different 
methodologies to produce identical results. The second 
thing is that for calculations other than our own in Ta- 
ble [Til different lattice constants or internal coordinates 
may have been employed, although it is our experience 
that, within reasonable limits, these never drastically re- 
duce the eccentricities (even though the actual numbers 
change) . 

In conclusion, we show that while the magnetic ground 
state of LaFeAsOi-^F^ can be obtained using any of 
several methodologies (all-electron, pseudopotential with 
ultra-soft or pseudopotential with PAW), the details of 
the magnetic properties, such as site magnetization and 
magnetic energy differences may be affected by pscudiza- 
tion of the crystal potential. More importantly, the 
results depend drastically on the exchange-correlation 
functional used (GGA vs. LDA) and on the position of 
As. The best results in terms of explaining the observed 
magnetic phase diagram are obtained within LDA (which 
is the recommended functional for itinerant magnets) and 
using the theoretical As positions. The fermiology and 
band structure are less sensitive to the details of the cal- 
culation, with the three all-electron codes we applied giv- 
ing identical results. Small Fermi surface details, such as 
the eccentricity of the central ellipses, differed between 
all-electron and pseudopotential codes and between dif- 
ferent pseudopotential codes, though a careful applica- 
tion of the latter yielded results that differ only slightly 
from the all-electron results. We have reported a refer- 
ence set of calculations for the basic electronic structure 
properties of the parent compound and for several dop- 
ings (both hole and electron). We believe these will pro- 
vide a consistent and accurate basis of knowledge upon 
which models and theories about this interesting com- 
pound can be built. 
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Note added in proof: After submission of this 
manuscript, we became aware of a paper by Ishibashi et. 
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a&. These authors applied an in-house pseudopotential 
code to four different fully AFM patterns (n > 1.8 fis 
on all sites) . Their stabilization energies in all four cases 
agree with our all-electron results within 7 meV/formula 
unit. Curiously, if one discards the FM pattern, the three 



energy differences among the remaining AFM states can 
be fitted with a two-parameter Heisenberg model fairly 
well. However, the fact that the FM ordering does not fit 
into this picture at all suggests that the fit is fortuitous. 
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